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Abstract. Spectral formation in steady state, spherical 
accretion onto neutron stars and black holes is examined 
by solving numerically and analytically the equation of ra- 
diative transfer. The photons escape diffusively and their 
energy gains come from their scattering off thermal elec- 
trons in the converging flow of the accreting gas. We show 
that the bulk motion of the flow is more efficient in up- 
scattering photons than thermal Comptonization in the 
range of non-relativistic electron temperatures. The spec- 
trum observed at infinity is a power law with an exponen- 
tial turnover at energies of order the electron rest mass. 
Especially in the case of accretion into a black hole, the 
spectral energy power-law index is distributed around 1.5. 
Because bulk motion near the horizon (1-5 Schwarzschild 
radii) is most likely a necessary characteristic of accretion 
into a black hole, we claim that observations of an ex- 
tended power law up to about m e c 2 , formed as a result of 
bulk motion Comptonization, is a real observational evi- 
dence for the existence of an underlying black hole. 

Key words: accretion — black hole physics — radiation 
mechanisms: Compton and inverse Compton — radiative 
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1. Introduction 

Spherical accretion into black holes as a Compton upscat- 
tering problem of low frequency photons was studied by 
Blandford and Payne (Blandford & Payne 1981a, hence- 
forth BP; Payne & Blandford 1981, henceforth PB). Espe- 
cially in PB the above authors solved the photon transfer 
equation in the case of a steady state spherically symmet- 
ric supercritical accretion into a central black hole with 
the assumptions of power law flow velocity and cold elec- 
trons T e = 0. An extension of the work of PB on accreting 
black holes to include the thermal motions of the elec- 
trons along with the bulk motion was performed by Colpi 



(1988). However, in the papers mentioned above, the au- 
thors solved the photon transfer equation under the as- 
sumption that the inflow extends up to the center (r = 0) 
and thus its thickness is infinite. 

Lyubarskij & Sunyaev (1982) solved the radiation 
transfer equation in the case of a finite optically thick 
medium. In a extension of the work of Blandford & Payne 
(1981b), they studied Comptonization in a radiation- 
dominated shock, taking into account not only the bulk 
motion of the electrons but also their thermal motion. 

In this paper we present the emergent spectrum from 
a spherical inflow by taking into account the effects of 
bulk motion and thermal Comptonization in addition to 
the boundary conditions of the inflow. For proofs and de- 
tailed discussions the reader is referred to our long pa- 
per (Titarchuk, Mastichiadis and Kylafis 1996, henceforth 
TMK). By extending the relevant work of Mastichiadis 
& Kylafis (1992, henceforth MK92) we allowed the inner 
boundary to be anything from fully reflective to fully ab- 
sorptive assigning the former case to a neutron star and 
the latter to a black hole. Also, we determined the eigen- 
values of the problem as functions of the mass accretion 
rate and we found numerically and analytically the general 
solution for the spectrum as well as its asymptotic forms. 
The obtained results in the black hole case show a remark- 
able similarity with the X-ray spectra of the Galactic black 
hole candidates in their high states. In §2 we present the 
solution of the radiative transfer problem while in §3 we 
give some examples and draw our conclusions. 

2. Radiative Transfer 

Consider a spherically symmetric accretion onto a com- 
pact object with mass rate M, where the bulk velocity of 
the infalling plasma is given by 



v(r) 



(1 - i)r s 



1/2 



(1) 
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with c the speed of light, r s the Schwarzschild radius and 
i the ratio of the escaping luminosity to the Eddington 
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one. The Thomson optical depth of the flow is defined by 



drn e (r)<jT =rn 



r{\-£) 



1/2 



(2) 



where n e (r) is the electron number density, <jt is the 
Thomson cross section, tti—M/Me and Me is the Ed- 
dington accretion rate defined by 



Me= 



L e AnGMnir, 



(3) 



Here Le is the Eddington luminosity, M is the mass of 
the central object, m p is the proton mass and G is the 
gravitational constant. 

Consider next low energy photons that find themselves 
in the accretion flow. As the photons diffuse outward, they 
scatter off the inflowing electrons gaining energy on av- 
erage. At some characteristic radius, called the trapping 
radius (Rees 1978; Begelman 1979), the conditions of op- 
tical depth and velocity of the inflowing electrons are such 
that the photons are advected inwards at the same rate as 
the one with which they diffuse outwards. Photons that 
either are emitted inside the trapping radius or find them- 
selves there cannot escape easily, but if they do, they gain 
significant amount of energy. 

Let us investigate next the relative importance of the 
bulk (converging flow) and the thermal motion of the 
electrons to the mean photon energy change per scat- 
tering. The mean energy gain per scattering of a pho- 
ton by thermal Comptonization < AE tc > is propor- 
tional to (v/c) 2 , i.e., < AE tc >« E(4kT e - E)/m e c 2 . 
On the other hand, the mean energy gain < AE c f > 
in the presence of a converging flow is proportional to 
v/c (e.g., BP). By using equations (1) and (2) we write 

< AE cf >« AEd{v / c) / d,T T = 4(1 - l)/m. Thus, we have 

< AE tc >/< AE cf > « (kT e - E/4)m/(l-£)m e c 2 < 
1/(5, where 



(kT e - E/A) m/(l 
e c 2 (1 - £) 



(4) 



This can be further written as 5 = 51.1 x T^ 1 m , with 
T 10 = fcT e /(10 keV). Hence, the bulk motion Comptoniza- 
tion dominates the thermal one if m T\q < 51. 

The dominant power-law and the high-energy cutoff 
of the spectrum can be understood in the following way: 
As it was shown above, the fractional increase in energy 
of a low-energy photon in its collision with accreting elec- 
trons of bulk velocity v(r) and temperature T e is given by 
< AE > lncr /E w 4(1 -£)/m + AkT e /m e c 2 . At the same 
time, the recoil effects cause a fractional decrease in the en- 
ergy of the photon given by < AE >decr/E w —E/m e c 2 . 
When E <C m e c 2 , the recoil effect is negligible, result- 
ing in a pure power-law spectrum. At high energies the 
two effects are comparable and the turnover in the spec- 
trum occurs at E c /m e c 2 w 4(1 — l)/ih + AkT e /m e c 2 , or 
E c /m e c 2 « 4(1 - £) rrT 1 when kT e < m e c 2 . 



The radiation spectrum observed at infinity can be 
found by solving the equation of radiative transfer for the 
photon occupation number n(r, v) (Eqn. [18] of BP). We 
remind the reader that the BP equation docs not take 
fully into account general relativistic effects; also the as- 
sociated diffusion coefficient has been derived in a fully 
self consistent manner. Both of these issues have been re- 
cently discussed (Zannias & Titarchuk 1996, henceforth 
ZT96). 

The two boundary conditions of the radiative trans- 
fer equation can be formulated in terms of the spectral 
energy flux F (cf. Eqn. [21] of BP). The first is that the 
total spectral flux integrated over the photosphere should 
depend only on E as tt — > or r — » oo. The second bound- 
ary condition is that we have a boundary with albedo A at 
some radius r b . The net energy flux through this surface 
is 



F(r b ,x) = -x 3 



I- A 

YTa 



(5) 



where x = hv/kT e and T e is the electron temperature 
considered to be constant throughout the flow. From the 
setup of the problem, r b is the smallest value that the 
variable r can obtain. 

If the inner boundary is fully reflective (A = 1), all 
input photons escape. This is the problem considered in 
MK92 who pointed out the analogy of the fully reflective 
boundary to a neutron star surface. On the other hand, if 
the inner boundary is fully absorptive (^4 = 0), then we 
have the equivalent of a black-hole horizon. However, since 
we do not take an infinite atmosphere but we consider in- 
stead a finite flow extending up to r b , our case is different 
from the one considered in PB and Colpi (1988). Rela- 
tivistic effects close to the horizon of the compact object 
have been ignored. 

The radiative transfer equation, satisfying the above 
formulated boundary conditions is solved in TMK. This 
solution is extended in ZT96 by taking into account rel- 
ativistic corrections. The spectra are calculated numeri- 
cally by using the iteration method to solve the boundary 
problem of the appropriate Fokker-Planck equation. The 
analytic approximation found by the separation of vari- 
ables method describes satisfactorily the emergent spectra 
in the whole energy range up to the high energy turnover 
E c ~ m e c 2 / m. Alternatively, as it is shown by TMK, 
this approximation can be written as a sum of two com- 
ponents: The low frequency injection is presented by a 
coherent solution and the high energy tail is presented by 
the fundamental mode. This mode can be presented by a 
convolution of the low-frequency source spectrum and the 
Green function I(x,xq) given by 



I = 



2n.x \x 



a+3+6 



(6a) 
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for x < xq and by 
1= b 



2fix T(2fi) 



Go) So 



-*(a; + t) a+3+tf t Q - 1 dt 



for X > Xq. 

Here b = a(a + 3 + 5), fx = |[(<5 - 3) 2 + 4 7 ] 



(66) 



1/2 
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2A (5, A is the first eigenvalue of the appropriate boundary 
problem and the spectral index is determined by 



a = fj, - (3 + 6)/2. 



(7) 



The normalization of I(x,x ) is chosen in such a way as 
to keep the photon number equal to l/x . 

As we noted earlier, the parameter which shows the 
importance of the bulk Comptonization effects relative to 
thermal Comptonization is S, given by Eqn. (4). In the 
case where kT e -C m e c 2 and to~ 1, then <5 > 1 and this 
leads immediately to an asymptotic relation for the energy 
spectral index a 

a = 2A 2 - 3 . (8) 

As long as S 3> 1, bulk Comptonization effects dominate 
the thermal ones and the value of a is independent of the 
temperature of the electrons. 

This boundary condition, along with the radiative 
equation implies that the eigenvalues are the roots of the 
equation (TMK) 

5 e \ */ ,2 5 7 . 5-2A 2 7 9 

2 " 3 Tb ) *<" A + 2 ' 2 ' n) + ^^rM-^ 2 , t>) = 

(9) 

(10) 



where = 3 m /2 and 



e = a + 3 — 



3(1 -A) 



2(1 + A) 



n 



r s (l-i) 



1/2 



By using Eqns. (9-10) it can be demonstrated (TMK) that 
the spectrum becomes harder as A increases. The depen- 
dence of the spectral index on the albedo A is displayed in 
Figure 1. We point out that for A — the power-law index 
a is always larger than 1 and approaches 1.5 for to ^> 1. 
For the fully reflective inner boundary case (A — 1), a — ► 
for to 3> 1. Thus, the hardest spectral slope (a = 0) cor- 
responds to the fully reflective inner boundary as it was 
shown in MK92. 

Of special interest is the case of a fully absorptive 
boundary A = 0. The obtained in this case spectral en- 
ergy power-law index a is shown in Figure 2 as a func- 
tion of to. As it can be seen, the spectral index is a weak 
function of the electron temperature in the wide range of 
mass accretion rates rh= 1 — 20 and electron temperatures 
kT e = — 3 keV. The temperature dependence of the spec- 
tral index becomes stronger for very high accretion rates 
to> 20 and for electron temperatures kT e > 10 keV. 

A second important consequence of the solution is 
that, despite the presence of the exponential in front of 
the integral in equation (6b), the power law can extend 



MASS ACCRETION RATE 



Fig. 1. Plot of the energy spectral index a versus m for dif- 
ferent albedo values A. Here kT e — 1 keV, i = and A = 
(solid line), A — 0.5 (dashed line), A = 1 (dash-dotted line). 
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Fig. 2. Plot of the energy spectral index a versus rh for dif- 
ferent electron temperature values kT e . Here A = 0, £ = and 
kT e = keV (solid line), kT e =0.1 keV (dashed line), kT e = 1 
keV (dash-dotted line) and kT e = 3 keV (dotted line). 

to x > 1 (see TMK). Thus, for kT e < to 6 c 2 , the power 
law extends to energies of order m e c 2 independently of the 
electron temperature. Only for energies x > S the spec- 
trum exhibits an exponential turnover. 

Figure 3 shows a comparison between the analytical 
approximation where the hard tail of the specrum is pre- 
sented by Eqn. (6b) and the numerical solution. The an- 
alytical approximation represents the emergent spectrum 
quite reliably in the whole energy range below the ex- 
ponential turnover. In this region two effects, the bulk 
motion upscattering and the Compton (recoil) downscat- 
tering, compete forming the hard tail. Because the ana- 
lytical representation does not take into account properly 
the inner boundary condition at energies E > E c there 
is a significant deviation of the analytical curve from the 
numerical solution close to the cut-off. When one takes 
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Fig. 3. Plot of the emergent spectral energy flux versus pho- 
ton energy calculated by numerical and analytical methods for 
kT e = 2.5 keV, m= 2, £ = 0, A = and kT bb = 0.833 keV. The 
solid line represents the numerically calculated spectrum while 
the dashed line represents the analytic approximation [the hard 
tail of the spectrum E > 15 keV is described by Equation (6b)]. 
The space source distribution is given by equation (11) with 
5*0 = 1 and /i = 0.3. 

into consideration the relativistic corrections, the hard tail 
becomes flatter, since Compton downscattering is less ef- 
ficient at high energies (ZT96). 

In Figure 4 the spectra are presented for different mass 
accretion rates. Here (as in Fig. 3) we have assumed exter- 
nal illumination of the converging flow by the low-energy 
black body radiation of an accretion disk having a char- 
acteristic temperature X^; furthermore we have assumed 
a space source distribution of the form (see, for example, 
Sobolev 1975) 

S(r) =5 r- 2 exp(-r T (r)/ Ai ) , (11) 

where So is a normalization constant and \i is the cosine 
of the angle of incidence at r ~ r s . All spectra are pure 
power laws in the high energy range from 15 kev up to 
the exponential turnover, which occurs at energy E c « 
(1 — £)m e c 2 / m. It is worth pointing out that, despite the 
fact that the accretion rate is taken to vary from 2 to 
10 Me, the resulting energy spectra have indices around 
1.5. 

3. Conclusions 

We have solved numerically and analytically the radia- 
tive transfer equation in the diffusion approximation in 
the case of thermal electrons infalling spherically onto a 
compact object. We have considered a finite medium (as 
opposed to an infinite medium considered in the related 
work of PB and Colpi [1988]) taking as inner boundary a 
totally absorptive or reflective surface, in the case where 
the compact object is a black hole or a neutron star, re- 
spectively. We have shown that the escaping spectral en- 
ergy flux exhibits an extended power law at high energies 




ENERGY (keV) 

Fig. 4. Analytical spectra for different mass accretion rates 
and kT e = 2.5 keV, t = 0, A = 0, kT bb = 0.833 keV. The m 
values are 2 (solid line), 3 (dashed line), 5 (dot-dashed line) 
and 10 (dotted line). The space source distribution is given by 
equation (11) with So = 1 and fj, = 0.3. 

up to a few hundred keV and has an index that depends on 
the inner boundary condition. The important parameter 
(in addition to the albedo introduced through the inner 
boundary condition) which determines the spectral slope 
is S (eq. 4) and as long as this remains much greater than 
unity the upscattering is mainly due to the bulk motion 
of the electrons. Therefore we find that, for a wide range 
of values of the dimensionless accretion rate m and of the 
electron temperature T e , the spectral slope is insensitive 
to both of these quantities. 

The above conclusion is especially true in the case of 
accretion into a black hole (absorptive boundary). The 
dominant energy spectral index takes values around 1.5 
for a wide variety of ra ranging from 1.5 to 10 (see Figs. 
2-4). Furthermore, this spectral index appears to be in- 
dependent of the electron temperature as long as this re- 
mains below 10 keV (see Fig. 2). Therefore we propose that 
our calculations may offer an explanation for the spectra 
observed by black-hole candidates in their high states (e.g. 
Ebisawa et al. 1993, Sunyaev et al. 1994 - for a recent re- 
view see Tanaka & Lewin 1995). The situation is changed 
only when either the temperature of the electrons is higher 
than 10 keV or when the escaping luminosity approaches 
the Eddington limit. In this last case the parameter 6 is re- 
duced due to the radiative force. In both cases the spectra 
are dominated by thermal Comptonization and are related 
to the hard spectra of black holes and neutron stars. 
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